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Charge separation in light-harvesting complexes occurs in a pair of tightly coupled chlorophylls at the heart of 
photosynthetic reaction centers of both plants and bacteria. Recently it has been shown that quantum coherence 
can, in principle, enhance the efficiency of a solar cell, working like a quantum heat engine (QHE). Here, we 
propose a biological quantum heat engine (BQHE) motivated by Photosystem II reaction center (PSII RC) to 
describe the charge separation. Our model mainly considers two charge-separation pathways more than that in 
the published literature. The two pathways can interfere via cross-couplings and work together to enhance the 
charge-separation yields. We explore how these cross-couplings increase the current and voltage of the charge 
separation and discuss the advantages of multiple pathways in terms of current and power. The robustness of 
the BQHE against the charge recombination in natural PSII RC and dephasing induced by environments is also 
explored, and extension from two pathways to multiple pathways is made. These results suggest that nature- 
mimicking architectures with engineered multiple pathways for charge separations might be better for artificial 
solar energy devices. 

PACS numbers: 42.50.Gy, 42.50.Nn, 84.60.Jt, 82.39.Jn 


I. INTRODUCTION 

Photosynthesis begins with the absorption of a photon 
which creates an excited state on a pigment molecule. The 
excitation is transferred between the pigments of light¬ 
harvesting complexes until it arrives at a reaction center (RC) 
in the pigment-protein complex, where the photon energy is 
used for the later dark stages of charge separation and con¬ 
version of energy from physical into chemical one. The effi¬ 
ciency of the energy transfer is very high, with a near unity 
yield H. This sparks the long-standing and increasing in¬ 
terest in the understanding of the physics behind the energy 
conversion within photosynthesis 12]. Recently, much atten¬ 
tion has been paid to the role of quantum coherence: growing 
experimental evidence M and theoretical models il0t - l22l] 
reveal that quantum coherence contributes beneficially to the 
high efficiency. Understanding the underlying mechanism of 
such natural system can assist us in designing novel nanofab- 
ricated structures for quantum transport and optimized solar 
cells. 

Viewing the photosynthetic reaction center as a quantum 
heat engine (QHE), the authors of Ref. i^ analyzed the 
charge separation in light-harvesting complexes. This treat¬ 
ment bridges the two seemingly unrelated effects attributed 
to quantum coherence in natural(photosynthesis) and artifi- 
cial(photovoltaics) light-harvesting systems. The common 
ground between photovoltaics and photosynthesis has also 
been investigated recently in |2^ . In analogy with a 
continuous Carnot-like cycle, Dorfman et al. showed that the 
power of a photocell based on Photosystem II reaction cen¬ 
ter (PSII RC) can be increased by 27% attributed to noise- 
induced quantum coherence—^Fano interference ll^ 1^ . 
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which was found in artificial photocells and lasing without 
inversion. Creatore et al. proposed the other mechanism and 
showed that the dipole-dipole interaction between two neigh¬ 
boring electron donors play a key role in enhancing the cur¬ 
rent and power of the photocell ll2a] . They claimed that this 
increasing can be up to 35%. 

In these studies, only one charge-separation pathway is con¬ 
sidered. Identifying the primary electron donors and dominat¬ 
ing charge-separation pathways has been a question of recent 
extensive research and debate. At the moment, there is much 
evidence that two main pathwa}^ make significant contribu¬ 
tion under ambient condition ll29l - [^ . This motivates us to 
ask the question that whether a multi-pathway scheme is more 
beneficial to designing artificial light-harvesting devices? In 
this paper, we will answer this question via numerical simu¬ 
lation for the current and power as well as the robustness of 
the current and power against the charge recombination and 
dephasing. 

The remainder of the paper is organized as follows. In Sec. 

II, we introduce a model of a quantum heat engine inspired 
by PSII RC to describe the charge separation, which includes 
a second pathway in comparison with the model proposed in 

im. A master equation describing the evolution of the 
two-pathway system is also derived in this section. In Sec. 

III, the concepts of effective voltage V and power are intro¬ 
duced and used to characterize the charge separation, new re¬ 
sults due to the second pathway are discussed. In Sec. IV, 
we focus on how the multiple pathways affect the behavior 
of j - V and P — V features. We explain the advantage of 
multiple pathways over simply increasing the decay rates and 
cross-couplings. In Sec. V, we take the charge recombina¬ 
tion and dephasing that might coexist in natural PSII RC into 
consideration. Numerical results of one- and multi-pathway 
models are shown, from which we further see the advantage 
of multiple pathways. Sec. VI is devoted to concluding re¬ 
marks. 
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II. MODEL 

To present the model, we first illustrate the structure of the 
Photosystem II reaction center complex in Fig.[T] The six pig¬ 
ment molecules are closely spaced in particular positions and 
orientations, and coupled by the dipole-dipole interactions re¬ 
sulting in exciton states. These pigments are distributed in 
two branches of protein matrix: D\ and D2. Pd\ and Pd 2 , a 
special pair of coupled chlorophylls, is located at the center 
of the PSII RC. They contribute mostly to the lowest energy 
states and are the primary electron donors, forming two exci¬ 
ton states which are denoted as |fli) and |fl 2 )- These two pairs 
of molecules are also coupled to the accessory chlorophylls 
Accdi and Accm located in the two different branches D1 
and D2, respectively. Phcoi and Pheoi are the remaining two 
pheophytin pigments coupled to the rest of the molecules, and 
act as the electron acceptors. Photosynthetic antennae trans¬ 
fer energy of the absorbed solar photons to the RC where the 
transmembrane charge separation takes place. Charge separa¬ 
tion in the core of pigment-protein RC complexes is the first 
energy conversion step in photosynthesis. The consequent 
electrochemical potential drives a chain of chemical reactions, 
including the reduction of NADP to NADPH, the synthesis of 
ATP and the oxidized part of the RC splits water, releasing 
molecular oxygen, all these leading eventually to the stable 
storage of solar energy. 

The photosynthetic reaction center may be analyzed as a 
biological quantum heat engine(BQHE) that transforms high- 
energy thermal photon radiation into low-entropy electron 
flux. As shown in Fig. |2] (a), a five-level QHE scheme de¬ 
scribes the photoinduced charge separation between the donor 
and the acceptor molecules interacting with thermal light. 
State \b) represents the lowest energy configuration in which 
both donors and acceptor molecules are in the ground state. 
States |fli) and |fl 2 ) correspond to single-exciton states in 
the first and second donors, respectively. Icr) is the charge- 
separated state with the excited electron transferred from the 
donors to the acceptor molecules. And \/3} describes the ion¬ 
ized state in which the system is positively charged (the ex¬ 
cited electron is assumed to have been used to perform work). 
The heat engine cycle begins with the excited electron pro¬ 
moted from \b) to |fli) and/or \a 2 ) after the absorption of a 
solar photon. This excited electron is then transferred to the 
acceptor by emission of a phonon (transition from |ai) and/or 
\a 2 ) to \a), corresponding to the formation of a radical pair 
P^^Acc ^2 in Fig-Ill*- Furthermore, the electron can be released 
from state \a) with a rate F, resulting in a current from lo-) to 
I/I) driving a chain of chemical reactions, leading eventually 
to the stable storage of solar energy. The current j - eTpaa is 
thus determined by the relaxation rate F and the population of 
|a). The acceptor-to-donor charge recombination represented 
by the decay rate - ;^fF, with x a dimensionless frac¬ 

tion, brings the system back to the ground state \b) but does 
not prod uce current, limiting the power output of our QHE 
0^ 1^ . Finally, to complete the cycle, we allow another 
population transfer to take place, emitting a phonon with ex¬ 
cess energy, bringing the electron back to the neutral ground 
state \b) with a rate F^. We also assume that the donor ex¬ 


cited states \ai} and \a 2 } undergo a dephasing process (via a 
rate rjep), which will influence the current and power gen¬ 
erated. This model including one charge-separation pathway 
(In') —> \j3)) can exhibit noise-induced quantum coherence due 
to Fano interference, which originates from the coupling of 
two levels to a common reservoir and thus the effect can be 
revealed by the cross-couplings. Fano interference can mini¬ 
mize the acceptor-to-donor charge recombination by inducing 
coherence between |fli) and |a 2 )- This brings about enhanced 
current and power output as we will show in section II. 

Although electron transfer in the PSII RC has been thor¬ 
oughly investigated and several charge-separation pathways 
that involve the formation of several different radical pairs 
were identified, works on the effect of multiple charge- 
separation pathways are rarely found. At the moment, much 
evidence shows that two main pathwa ys m ake significant con¬ 
tributions under ambient conditions 0291 - 1^ . Therefore we 
emphatically study two pathways, i.e., we add another accep¬ 
tor molecule that consists of a charged-separated state and an 
ionized state to the model in Fig. |2(a). As shown in Fig. 
(b), the initial excitation of states |ai) and |a 2 ) can be trans¬ 
ferred to the two acceptor molecules in states lai) and |a 2 ), 
respectively, with the excess energy radiated as a phonon, and 
further produce electric current and returning back to \b) via 
^i) and/or ^ 2 )- Here, the total current in the two pathways is 
given by j - eFipaja, +eF 2 Pa 2 Q.,. The structures ofRCs from 
green plants, algae, and bacteria differ in the nature and pre¬ 
cise orientation of the constituent pigments, resulting in dif¬ 
ferent spectroscopic and dynamical parameters. Nevertheless, 
our six-level scheme may be applied to all RCs as discussed 
below. 

With this knowledge, the dynamics of our heat engine can 
be described by the following master equation, 

^--i[//5,Ps]+Tps, (1) 

at 

where Hs - 2 \j) 01 is the system Hamiltonian, 

j=b,a\,a2, 

Q'l,Q'2,y3l,j02 

with Ej the energy of the /th level. The superoperator L de¬ 
scribing the effect of reservoirs and acceptor-donor recombi¬ 
nations can be decomposed as 

Fps -LhPs + LclPs + Lc2Ps + LrelPs 
+ Ljepps + Lrecps ■ 

For the high temperature photon reservoir (HTR), its effects 
can be described by 

Ups = ^ ^[(nh-\-l)(o-kipso-l^ + o-hi,pso-l^ 
k,l=l,2 

t ' t \ / t 

- cr I,i^pg - PscrJ^^cr 1,1) + ni,{cr^^j^ps cr 

+ CrliPs^^hk - O'hkO'liPs - Ps^'hlO'li;)], 

where &hk = \b){aic\ (k - 1,2), and Hf, is the average number 
of solar photons, {k - 1) denotes the decay rate, where 
Jnh - 7\h and 722 ^ = T 2/1 represents the spontaneous decay 
from the upper level |ai) and |a 2 ) to level \b), respectively. 
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FIG. 1: (Color online) Arrangement of six core-pigments in the PSII RC. It consists of the special pair Pdi, Pd 2 , the two accessory chlorophylls 
Accdi, Accd 2 , and the two pheophytins Pheni, PheD 2 . 



(C) 


FIG. 2: (Color online) Schemes of the BQHE model based on the photosynthetic reaction center. HTR denotes the high-temperature photon 
reservoir, while LTRl and LTR2 stand for the low-temperature phonon reservoirs. HTR induces transition from the ground state \b) to the 
single-exciton states |fli) and/or \a 2 ). LTRl induces transition from |fli) and/or |fl 2 ) to the charge-separated state |a) (|ai) and |q' 2 ) for (h), In',-) 
(i = 1,2, 3...) for (c)). LTR2 induces transition from the ionized state \p) (|j 8 i) and 1 ^ 82 ) for (h), |j 8 ,) {i = 1,2,3...) for (c)) to the ground state. 
The three models differ in the number of donor, i.e., pathway. 
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The noise-induced quantum coherence is closely related to 
the terms with {k + 1), which we will refer to as cross¬ 
couplings. They describe the effect of interference, and we 
will assume yn/, = yiih, and juh = r]h ^/rki^yih (k + 1) where 
rjt, - I represents the maximal coherence and rjh — 0 the min¬ 
imal coherence. It is worth noting that rji, - 0 represents a 
case in which the two donor molecules interact with their own 
(independent) reservoirs. Thus there is no coherence induced 
by noise ll^ in this situation. 

Superoperator Lc 2 corresponding to the second low temper¬ 
ature phonon reservoir (LTR2) has the form 


Lrei describes a process that the system in state lori) (|£r 2 )) 
decays to state ^i)(|yS 2 ))- It leads to the electronic current 
proportional to the relaxation rates Ti and r 2 as defined later. 


p 

LreiPs = ^ -:^{\f^k}{ak\ps \a!k)<J^k\ 
k=l,2 ^ 

- \ak){ak\ps -ps \ak){ak\). 


(6) 


We describe the dephasing of the system on states |ai) and 
|fl 2 ) by 


LciPs = 2j + o-ckps 0 - 1.1 

k,l=l,2 

t t t ^ (4) 

- o-'^/O-ckPs - PsCTck^ci) + nc 2 {cr[kPsO-ci 

+ o-\ps O'ck - o-ckO-liPs - Ps O' c/crjj.)], 

where &ck = \b)(J3k\ {k - 1,2) denotes the lower operator. 
nc 2 is the average phonon number of the cold (phonon) reser¬ 
voir LTR2. Tkic (k - 1) is the corresponding decay rate, with 
Tile = Tic and r 22 c = r 2 c representing the spontaneous decay 
from the upper level |y6i) and |/32) to level \b), respectively; 
and Tkic (k F) is a cross-coupling that describes the effect of 
interference, with ri 2 e = r 2 ic. In numerical simulations, we 
will set r kic - TJcn VTicTic {k + 1) with rje^ - 1 representing the 
maximal coherence and tJc^ - 0 the minimal coherence. 

The quantum coherence induced by the low temperature 
phonon reservoir (LTRl) can be described by, 

U\PS = + l){OklPsoi„ 

kj,m,n=\,2 

"I" ^wiips ^i^i ~ ^j^f^mnPS ~ PS^\in^kl) (5) 

+ (T]cl + ffifi — 0~fu,i(T 

- PsOkiol,„)], 

where &ki = \ak){ai\ (k,l = 1,2) or o-,„„ = |a,„> (a„| (m,« = 
1,2), and iid is the average phonon numbers of LTRl. 
Jkimnikl - mn) are decay rates from the upper levels |ai) and 
|fl 2 ) to the lower levels \a\) and \a 2 ), respectively. For kl + mn, 
terms with ykimn would induce quantum coherence leading to 
the interference. Note that the cross-couplings here are com¬ 
plicated more than those of HTR or LTR2. There are two 
types of cross-couplings due LTRl. One is the same as that 
in previous works lUlUl, defined as ykimn = nl, iykiymn 
{kl + mn), including ymi and yi 222 - These two cross¬ 
couplings represented by terms with yi 2 e in the model shown 
in Fig.|2](a). The other one is defined as yktmn - rfX y/ykiymn 
{kl + mn), including terms with yiii 2 , yii 22 , ri 22 i, 72122 - 
These terms are not considered in the model depicted in Fig.|2 
(a). The first type of cross-couplings couple \a\) and \a 2 ) with 
the same lower lever. Namely, terms with yii 2 i couple |ai) and 
\a 2 ) to Icri), while terms with yi 222 to \a 2 )- In contrast, the sec¬ 
ond type of cross-couplings, the lower levels are different. In 
the later numerical simulations, we will choose - 1 repre¬ 
senting the maximal coherence and 0 the minimal coherence. 
We will discuss the effects caused by these cross-couplings in 
section IV. 


Z F ifgp 

-—{\ak){ak\ps \ak){ak\ 

k=l.2 

- \ak) {ak\ps -ps \ak}{ak\). 


(7) 


with Tgigp being the dephasing rate. 

Finally, the loss channel due to acceptor-to-donor charge 
recombination is described by. 


LreePS = Y lOk) (b\ 

k=l,2 

- \ak) {oklps - Ps IQ'a) (^aD- 


( 8 ) 


Here the dimensionless x stands for the recombination rate. 


III. EFFECTS OF CROSS-COUPLINGS 


TABLE I: Parameters used in the numerical simulations. 


Parameters 

Units 

Values 


cm' 

14,856 


cm' 

1,611 

r,h 

cm' 

0.01 

riin{=72c) 

cm' 

140 


cm' 

18 

r.(=rj 

cm' 

200 

r,(=r) 

cm' 

1,000 

E dep 

cm' 

41 



60,000 



4.4x10"'‘ 

Ng 


4.4x10"^ 

Z 


0.2 
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In the present system, the current generated can be thought 
to flow across a load connecting the acceptor levels |a) and \[i). 
Introducing an effective voltage V as a drop of the electrostatic 
potential across the load, we obtain eV - Ea—Ep+ksT In(^) 
for the model in Fig.|2](a), where e is the electric charge and 
Ei is the energy of state |/). 

In our model shown in Fig. |2] (b), choosing Ei^ - E,^ 
(i = a, a,13) and setting all parameters the same for the two 
pathways, we have eV = E^^ - £/?, + kgTlni^y), and 

eV = ln(^^^), which define the voltage V. The 

current, accordingly, is defined by j - eT\pa^al + eFap^a,. It 
is similar to the parallel circuit in classical electromagnetism. 
Based on the current and the voltage, we easily obtain the 
power output P - j-V. We apply this to the steady-state of the 
system to characterize the performance of our BQHE and the 
photovoltaic properties of the complex, i.e., the steady-state 
current-voltage (j - V) and power-voltage {P - V) character¬ 
istics. Using the steady-state solution of Eq. ([1]), we plot the 
j-V curve and power with different F, while keep the other 
parameters fixed. Note that F ^ 0 (j ^ 0) corresponds to the 
open-circuit case, and in the short-circuit case, V —> 0. 

The parameters used in our simulation are listed in Table 
I. These parameters are chosen from recent literature, e.g ., 
and they are used in the simulation in 
The energy differences are defined as Ea/b - Ea^ - Eg, Ejj.g - 
E/j. - Eh and Ea/aj = - Eaj, where i,J =1,2. Parameters 

in the brackets of Table I are used in the one-pathway model 
in Fig. 12 (a). In the following, we will discuss the effect of 
cross-couplings separately. 


A. Cross-coupling Tnc induced by LTR2 



Time[fs] 


Time[fs] 



10 _ 
8 ^ 
6 2 
4 ' 

2 
0 


< 1 , 


Voltage (V) 


Voltage (V) 


FIG. 3 : (Color online) Fig.[ 3 ](a)-(d) show the time evolution of the 
population on each level of two-pathway model corresponding to 
Fig.| 2 ](b) at room temperature T = 300 K: pj,;, in red-dashed, Pa^a^ 
(Pa2a2) la orange-dot-dashed, pa^ai ip 112112) in green-solid, and p^j^j 
(Pfiifii) in hlack-dotted line. The cross-coupling Fi2c = rjn VTuF^ 
takes 77^2 = 0 , 0 . 3 , 0 . 6,1 for (a), (h), (c) and (d), respectively. The 
other cross-couplings 7,2* = pi, VriTXM, ykinm = i 7 ki 7 mn ikl 
mn) and jkinm = 'icj ^fyuTnm (kl inn) take the maximal values, i.e., 
'7/1 = '7*1 - vii - )■ (®) nnd (f) are plotted for the current and power 
generated as a function of effective voltage V at room temperature: 
t]c2 = 0 in solid, 77^2 = 0.3 in dotted, 77^2 = 0.6 in dot-dashed and 
77^2 = 0.9 in dashed line. 


Interestingly, the j-V and P - V behaviors demonstrate 
that as cross-coupling strength F 12 c increases, the current and 
power decrease, see Fig. [2(e) and (f). This is different from 
the results in previous works Note that, the param¬ 

eters chosen for the two pathways are the same (see. Table 
I), thus the population on |ai) and |a 2 ) is the same, so is the 
population on lai) and |q' 2 ), ^ 1 ) and ^ 2 )- From Fig. |2(d), we 
see that, the maximal Fi 2 e (r 7 c 2 = 1) benefits the population 
on ^ 1 ) and ^ 2 ) and leads this population dominant over the 
others. Therefore we cannot plot j-V characteristics in this 
case as = Pa 2 a 2 - 0- It was shown in Ref. that 

the power is increased and the current-voltage characteristic 
of the heat engine gets better with Fano interference. Here 
Fano interference characterized by Fi 2 c, however, shows an 
opposite effect. We will explain this observation in subsection 
D. So, in the following numerical simulations, we will take 
ri2c = 0. 


B. Cross-couplings induced by LTRl 

Fig.|4|shows that, current and power are enhanced as 711,21 
and 712,22 increase, which coincides with the results predicted 


by Dorfman et al. This gives rise to a question that why does 
Fi 2 c exhibit effects different from that by 711,21 and 712 , 22 ? 

Fig. 12 shows that, population in each state seems no dif¬ 
ference for any values of the second type of cross-couplings 
Jkimn - '7c, -sjykiymn (kl P mn). Thus different from couplings 
with 71121 ( 71222 ) and 712 c, the second type of cross-couplings 
has marginal effects on the current and power. This is anal¬ 
ogous to quantum beats, which can occur in a V type atom, 
but can not in a A type atom based on the theory of quantum 
electrodynamics (OED) ll^ . as shown in Fig.|2 V-type atom 
in excited states decays via the emission of a photon with fre¬ 
quency vi or V 2 . Since both emissions lead to the same final 
state, we cannot determine along which path, vi or V 2 , the 
atom decays. This uncertainty leads to an interference be¬ 
tween vi and V 2 , exhibiting quantum beats. A A-type atom 
also decays via the emission of a photon with frequency vi 
or V 2 . However, an observation of atom at a long time af¬ 
ter the emission can tell us which channel (1 or 2) was taken 
(atom in |77i) or |^' 2 ))- Consequently, we expect no beats in 
this case. For the model of Fig. [2] (b), transitions described 
by the second type of cross-couplings lead to different lower 
states, which do not induce interference. Thus the second type 
of cross-couplings has no effect on the j-V and P-V behav- 
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FIG. 4 : (Color online) Fig.| 4 ](a) - (d) show the time evolution of the 
population on each level in the two-pathway model shown in Fig. 
(b) at room temperature T = 300 K. pi,i, was plotted in red-dashed, 
Poiai (Pa2m) in orange-dot-dashed, green-solid, and 

(P/S2/S2 ) in black-dotted line, in the cross-couplings ykimn = 
Pci ^yWYmn {kl t- mil) takes = 0 , 0 . 3 , 0 . 6,1 for (a), (b), (c) and 
(d), respectively. The other cross-couplings 712;, = p;, yjyihjih and 
yumn = Pc, ^llkiynm {kl 4 - mil) take the maximal values, i.e., rji, = 
77^1 = 1 . As mentioned in the last subsection, F|2c = 0 . (e) and 
(f) are plotted for the current and power generated as a function of 
effective voltage V at room temperature. Dashed line is for 77*^ = 0 , 
dot-dashed for 77'^ = 0 . 3 , dotted line for 77*^ = 0 . 6 , and 77'^ = 1 is 
ploted in thin-solid line. 


FIG. 5 : (Color online) Fig.[^(a)-(d) show the time evolution of the 
population on each level of the system depicted in Fig.| 2 ](b) at room 
temperature T = 300K: pbt, in red-dashed, pnjc, ipaiai) in orange- 
dot-dashed, p„,„, ip0,2^2) in green-solid, and p^,^j i” black- 

dotted line. The cross-couplings yumn = VwXmn {kl + mii) take 
= 0,0.3,0.6, 1 for (a), (b), (c) and (d), respectively. The other 
cross-couplings 7,2* = rjh iyihyih and yumn = pi, iyuymn {kl + mn) 
take the maximal values, i.e., 77* = Pi, = 1- r,2c = 0. (e) and (f) are 
plotted for the current and power as a function of effective voltage 
V at room temperature T = 300K. Actually, we plot four curves in 
(e) and (f) with 77^^ = 0,0.3,0.6, 1 , respectively, although all them 
overlap. Since (a)-(d) do not show any difference for different 77^^, 
the j — V and P - V behaviors are exactly the same for different 77^^. 


iors. 

With this consideration, in the following numerical simula¬ 
tions, we take 7i;„,„ = ^lykijmn = 0. 


C. Cross-couplings induced by FITR 

As expected, ynh benefits the transition as seen from Fig.| 7 ] 
Inspired by the physics behind the quantum beat, we claim 
that interference can play an important role in the current and 
power only when the upper two levels \a\) and \a2) are coher¬ 
ently excited at the initial time. When the coherence between 
\a\) and \a2) is absent, i.e., seting 712/, = 0, 712c itself may sup¬ 
press the exciton transfer. Equivalently, setting 712c = 0 , 712/, 
itself may inhibit the current and power. In terms of steady- 
state solution to Eq. [U we further confirm this observation, 
and it shows how these two types of interference work to¬ 
gether to affect the transitions. 


D. Combined effects of HTR and LTRl 


Eor simplicity, we plot the steady-state j - V character¬ 
istic and power for one-pathway model in Eig. 0 We set 
7i2c = Pc, ^jyxcyic and ynh - Ph yjyihyih, thus = 1 
{rjh = 1) corresponds to the maximal coherence and p'^ = 0 
{ilh - 0 ) the minimal coherence. We observe that, when yi2h is 
very small, the current and power decrease with the increasing 
of 7 i2c. Similarly, when 712c is very small, large 712;, inhibits 
the transitions. Recall the quantum beats, if one of ynh and 
7i2e is very small, |ai) and \a2) are almost independently ex¬ 
cited. Eor this reason, when one cross-coupling term becomes 
very large, the current and power will be suppressed. When 
both 7127, and 712c increase to a certain value, strong cross¬ 
couplings can enhance the current and power. This conclusion 
holds for 71121 and 71222 in the two-pathway model. Eor the 
terms with ri2c discussed in subsection A, ^1) and ^2) obvi¬ 
ously, are not in a coherent state. As a result, the transition 
process is suppressed as ri2e increases. 


(A 3 eOI) "'ll d 
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a. 


(a) 

V - type atom 


\ / 
\ / 

\ / 


b' 


\Wv(t))= Z ^i\b^)^i\b,\) + C2\b,\) 

i=ai ,02 ,b 


a 


(b) 

t \ 

/ \ A - type atom 

Vj ' ^2 





/=t7| ,02 


FIG. 6: (Color online) Three-level atomic structures for (a) V-type 
and (b) A-type quantum beat. 


IV. EFFECT OF MULTIPLE PATHWAYS 


Electron transfer in the PSII RC has been thoroughly stud¬ 
ied and several charge-separation pathways were identified. In 
this section, we investigate the role of multiple pathways on 
the behavior of j - V and P - V. The simulation results are 
presented in Fig.|9] From Fig.|9]we find that as the number of 
pathways increases, the BQHE produces stronger current and 
power. Especially, the model of two pathways, corresponding 
to Fig. |2] (b), shows a current enhancement of 76.8%, com¬ 
paring with the case of only one pathway shown in Fig.|2(a). 
Nevertheless, with the number of pathways increase further, 
the enhancement rate decreases. The current of three-pathway 
model shows 25.9% more than that of two-pathway model, 
and the current of four-pathway model just 16.7% more than 
that of three-pathway one. The observations can be under¬ 
stood as follows: the current and power of the BQHE sharply 
depends on the rate of transition from the ground state \b) to 
state \a\) or \a 2 ), i.e., the decay rate in HTR. More pathways 
mathematically equal to increasing the strength of the decay 
rate in FTRl and LTR2. When y\h, yih and ynh take a fixed 
value, the energy flux density is fixed. The strength of the de¬ 
cay rate in FTRl and LTR2, or rather the number of pathways 
should match the transport rate in HTR, otherwise, excessive 
pathways can not significantly improve the transition process. 
We further propose that in order to accelerate the transfer pro¬ 
cess, multiple pathways in the BQHE (or similar structure in 




FIG. 7 : (Color online) Fig. | 7 ] (a)-(d) show the time evolution of 
the population on each level of the system in Fig. (b) at room 
temperature T = 300 K. pj* is shown in red-dashed, Paiai iPa^ai^ 
in orange-dot-dashed, Pa^al (Paj^) in green-solid, and p^,^j (P/Jj/Jj) 
in black-dotted line. The cross-couplings 712;, = 77;, yjy\hy2h take 
Tji, = 0 , 0 . 3 , 0 . 6 , 1 for (a), (b), (c) and (d), respectively. The other 
cross-couplings 7*/,„„ = 77*^ s/ykiynm (kl + inn) takes the maximal 
value, i.e., 77 ^^ =1. y*/,,,,, = 77 ;, = Q (kl t mn) and Fi 2 c = 0. 

(e) and (f) are plotted for the current and power as a function of ef¬ 
fective voltage V at room temperature. 77* = 0 is plotted in dashed, 
77;, = 0.3 in dot-dashed, 77/, = 0.6 in dotted and 77/, = 1 in thin-solid 
line. 



FIG. 8: (Color online) Counter plot of the current as a function 
of cross-couplings 712c and 712*. As 712c = 77'^^ ^yicyic (ynh = 
r]h yjyihyih), v'ci = 1 (t/h = 1) corresponds to the maximal coherence 
and 77'^! = 0 (77h = 0) the minimal coherence. 
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FIG. 9 : (Color online) The current and power as a function of effec¬ 
tive voltage V at room temperature for different numbers of pathway. 
The results for one pathway are shown in dashed, two pathways in 
dot-dashed, three in dotted, four in thin-solid line. Cross-couplings 
corresponding to LTR 2 and the second type of cross-couplings cor¬ 
responding to LTRl are set to 0 . All the other cross-couplings take 
their maximal values. 


the natural photosynthetic reaction center) are a better choice 
comparing with simply increasing the strength of the decay 
rates and cross-couplings related to LTRl and LTR2 within 
one-pathway model. As shown in Appendix, the decay rates 

or cross-couplings are given by 'yi;/nm - . from 

which we know that these rates are determined by the ambi¬ 
ent phonon reservoir and the diploe moment of |fli, 2 ) \b} 

for LTRl (|y 6 i, 2 ) <-> 1^) for LTR2). The natural condition 
of reservoir is hot and wet, and generally unstable. And the 
structure of RC complexes is born fixed, resulting in constant 
diploe moment between different states of pigment molecules. 
Therefore, these factors cannot be utilized to enhance the re¬ 
action efficiency. Nevertheless, more charge-separation path¬ 
ways means increasing reaction channels. This seems easy 
to achieve. In the next section, considering acceptor-to-donor 
charge recombination and noise-induced dephasing, we will 
show that multiple pathways is still a better choice in contrast 
with increasing the strength of the decay rate in LTRl and 
LTR2. 


V. EFFECTS OF RECOMBINATION AND DEPHASING 

In this section, we examine how charge recombination at 
the acceptor/donor interface and dephasing of |ai) and \a 2 ) 
affects the current and power. In Fig. [10] (a), we plot the 
relative change of current as a function of the recombina¬ 
tion rate x- The relative change of current is defined by 
(jCt = 0 ) “ j(x))lj(x - 0 )- We calculate this relative change 
for both cases of maximal coherence and minimal coherence 
with one and two pathways, respectively. To determine which 
one, multiple pathways or large decay rates in LTRl and 
LTR2, is a better choice, we double these decay rates and 
the relaxation rate in one-pathway model with the other pa¬ 
rameters fixed. Fig. [T 0 |(b) shows the relative current change 
for both cases of maximal coherence and minimal coherence 
with one and three pathways, respectively. Similarly, we triple 
these decay rates and the relaxation rate in one-pathway model 
with the other parameters fixed. Fig. [T 0 |(a) and Fig. [TOl(b) 
show that as the recombination rate increases, the current 




FIG. 10 : (Color online) (a) is plotted for the current ratio defined by 
(j(X = 0 ) - j(x))l jix = 0 ) as a function of the recombination rate for 
one- and two-pathway models. The blue-solid and red-dashed lines 
are for the maximal coherence, i.e., all the cross-couplings take the 
maximal values, except that cross-couplings induced by LTR 2 and 
the second type of cross-couplings induced by LTRl are set to 0 . The 
blue-dot-dashed and red-dotted lines are for the minimal coherence, 
i.e., all the cross-couplings take the minimal values 0 . The blue lines 
stand for one-pathway model with the strength of the decay rates in 
LTRl and LTR 2 and the relaxation rate doubled and the other param¬ 
eters fixed. The red lines stand for two-pathway one. (b) is the case 
of one- and three-pathway models. The blue-solid and red-dashed 
lines are for the maximal coherence while the blue-dot-dashed and 
red-dotted lines are for the minimal coherence. The blue lines stand 
for one-pathway model with the strength of the decay rates in LTRl 
and LTR 2 and the relaxation rate tripled and the other parameters 
fixed. The red lines stand for three-pathway one. 


decreases monotonically as expected. Nevertheless, the set¬ 
up with maximal coherence can better reduce the impact of 
acceptor-to-donor charge recombination. Therefore, noise- 
induced coherence between |ai) and |a 2 ) due to Fano interfer¬ 
ence can diminish the influence of recombination. Besides, 
we observe that the larger number of the pathways is, the 
more robust the system against the recombination. There¬ 
fore, in the presence of the additional loss mechanism due to 
acceptor-to-donor charge recombination, the maximal coher¬ 
ence scheme/multi-pathway set-up is more efficient than the 
minimal coherence scheme/one-pathway one. 

Next we explore whether, in the presence of dephasing, the 
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multi-pathway set-up can still benefit the current and power. 
We calculate the relative change of current as a function of the 
dephasing rate Ydep and present the results in Fig.[TT] The rel¬ 
ative change of current has also been calculated for different 
numbers of pathways. Similarly, to compare the two choices 
of multiple pathways and increasing the strength of the de¬ 
cay rate in LTRl and LTR2, we modify these decay rates and 
the relaxation rate in one-pathway model with the other pa¬ 
rameters fixed, just as in the case of the relaxation rate dis¬ 
cussed above. Observing the results, we find that the current 
decreases monotonically as the dephasing rate increases. Be¬ 
sides, the more pathways, the more robust against the dephas¬ 
ing. Therefore, the multiple-pathway model is more efficient 
than the one-pathway one even in the presence of depasing 
process mechanism. 
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VI. CONCLUSION 

In conclusion, we have studied the dynamics of two¬ 
pathway biological quantum heat engine and calculated the 
steady-state current-voltage (j-V) and power-voltage (P- V) 
characteristics. We explored the effect of various cross¬ 
couplings, and found that only the cross-couplings describ¬ 
ing transitions between high levels and the same lower energy 
level have effects on the j — V and P - V behaviors. We ex¬ 
plained the physics behind this observation using the concept 
of quantum beat. Furthermore we find that the current and 
power can be increased due to Fano interference when the up¬ 
per two levels are in a coherent supposition. Finally, we calcu¬ 
late the j-V and P -V characteristics subjected to acceptor- 
to-donor charge recombination and dephasing, and show that 
noise-induced quantum coherence helps in the suppression of 
the influence of the acceptor-to-donor charge recombination. 
These results suggest that multi-pathway BQHE can benefit 
the charge separation and the light-harvesting yields, and it 
is a better choice in contrast with increasing the strength of 
decay rates and cross-couplings. 
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Appendix A: The derivation of Eq. jlj 

In this part, we will present a derivation for the master equa¬ 
tion used in the maintext. We consider a total Hamiltonian 

H = Hs+Hr+Hsr, (A1) 

where the system Hamiltonian 

Hs - htOh \b} {b\ + hoja, |fli> (ail + htOa, Wi) {ail, (A2) 

with htOa b - Ea b {i - 2) being the energy of states |a,), 

\b}.^ 

Hr represents the Hamiltonian of the high temperature pho¬ 
ton reservoir (HTR) 




FIG. 11: (Color online) (a) is the relative current change (defined by 
(j(Ydep = 0)- i{Tdep))ljiVdep = 0)) as a function of dephasing rate. 
The blue-dashed line is for one-pathway model with the strength of 
the decay rates in LTRl and LTR2 and the relaxation rate doubled 
and the other parameters fixed. The red-solid line is for two-pathway 
one. (b) is the case of one- and three-pathway model. The blue- 
dashed line is for one-pathway model with the strength of the decay 
rates in LTRl and LTR2 and the relaxation rate tripled and the other 
parameters fixed. The red-solid line is for three-pathway one. The 
dimensionless parameter^ = 0.2 as in Table I. 


= Z (A3) 

P 

Here and bp denote the creation and destruction operators. 
The interaction Hamiltonian in the rotating-wave approxima¬ 
tion (RWA) is 

Hsr = h'^gyp&hibl + h'^g2qd-h2bl + H.c., 

P 


with &h\ - \b){a\\ and &h 2 - \b){a 2 \ being the lower op¬ 


erators. 

Pa-rb'^q*^p 


The coupling constant gip = 
where pa^b {pa 2 b) is the dipole moment of 1^3 1 ) 




h 
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1^) (l« 2 ) 1^) )■ £/! (£q) is the polarization of the field. The Substituting V{t) into Eq. (IA6b . we have, 

electric field per photon is ■ 

It is convenient to work in the interaction picture. The 
Hamiltonian, in the interaction picture, is given by 


k 

+ + H.C.. 

q 


(A5) 


The equation of motion for the system density operator ps is 


dpsit) 

dt 


- ^Tr«[y(f),P5(0)®pR(0)] 
a 

- [V(t),[V(t'),ps(t')»PRm]dt'. 


(A6) 

I 


dps(t) 

dt 


- dt'\Y^gipgip, 


[atibp, [bl,&h2,Ps(t') ®P^(0)]] 

p,q' 

^ [d-h2bq, [^],,o-/,i,ps(f') ® P«(0)]] 

q^p' 


+ ^ 8lp82q’ 

p^q' 


+ 

q,p' 

+ ^ 82q82q' 
‘M’ 


,'^'(‘^“2-'',[<T^2fe„ [^^,<^ft2,Ps(f') ®P«(0)]] 

X [bl&Hu [&lbp,,ps{t') ®p«(0)]] 

P^P' 

2 [bl&hu [&l2bq',Ps{t') ®Pfi(0)]] 

[(T,tiV,ps(f') ®Pr(0)]] 


+ 

p^q 

+ ^ 82q8\p‘ 

qqy 


■ J] 82 q 82 q'e ''‘f'>‘'TTR[bl&h2,[&l2K'’PsiO®PR(0)]]\ 

q^q' ) 


(A7) 


The sum over p may be replaced by an integral through 


Z 



(A8) 





(A9) 


where Vpi,t is the volume. Neglecting all memory effects and 
assuming that the density matrix is a slowly varying function 
of time, i.e., p(t') ^ p(t), we obtain the integration over time 
as 


Substitute dASb and ( IA9b into ( IA7b and note that (hp'j - 
(hi) = 0, {bpbp,) = {blbl,) = 0, (^pV) = bpSpp', and 
{bpbl,) = {rip + 1)6 pp'. The master equation reduces to 
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i)(2o-mPs(0o-1i - o-IiO-MPs(0 -Ps(Oo-IiO-h)] 

+ P^slp [^P {^^liPs(t)^h\ - d-hi&l^psit) - psiOd-hi&li)] 

+ P^slp [(”p + 1)(2 o-/,2Ps(0<3'|,2 - ^l2^h2Ps{t) - Ps{t)^l2^h2)\ 

+ [”/> {'^^laPsiOd-hi - d-h2^l2Ps{t) - Ps{t)d-h20-\^\ 

+ P^glpg2p [{tip + l)(<TMPs(0o'l2 + O'A2Ps(0o'Ii)] 

+ P^glpg 2 p [ftp (o-|iPs(0o-/,2 + d-l 2 Psit)d-hi)] 

+ P^gipg 2 p [{rip + i)((T/,2Ps(0o'1i + 

+ P^g\pg 2 p [rip [&'l 2 Psif)^h\ + ^liPs(t)d-h 2 )]} ■ 


The indicator p here is redefined as p = Introducing notations 


2p^Vph,gkpgip 

Jklh - - 

nc 


(All) 


where k,l - 1 , 2 , gkjp is the coupling constant corresponding to transition with frequency Wa,,*, and the average occupation 
numbers of photons tip - hp, we rewrite Eq. (lAlOb in a simple form as 


dpsjt) 

dt 


kj=l,2 


+ cr hkPs O'!; “ cr/,^ps - Ps O' 1,^0-hi) 


■ rihicrlf^psCThi + crliPsCThk - o-hkO-\iPs -psO-/,/cr^j,)]. 


(A12) 


The cross-couplings ykih for k + I corresponds to the strength 
of the Fano interference, and jkkh - Jkh denotes the sponta¬ 
neous decay rate. With this consideration, we finally obtain 
Eq. a. 

By the same procedure outlined for the transition \b) <-> 
\ai 2 ) induced by the high temperature photon reservoir 
(HTR), we easily obtain the Lindblad term describing the tran¬ 
sition process \b} <-» |y 6 i, 2 ) induced by the second low tem¬ 
perature reservoir (LTR2) in Eq. (IHi, where &ck = \b}(J3k\ 
(k - 1 , 2 ) is the lower operator. nc 2 is the average phonon 
numbers of LTR2. The spontaneous decay rates and the cross¬ 
couplings are connected by, 

p _ 2p Vphngkpglp (A13) 

nc 


where the indicator p is defined as p = gk,ip is the cou¬ 
pling constant for transition with frequency a>p^ Vpim is the 
phonon volume. 


For the transition \a\ 22 ) <-> pi, 2 ) induced by the first low 
temperature phonon reservoir (LTRl), derivation is tedious 
but similar to the presented one. The complicity comes from 
the extra charge-separated states. The Hamiltonian describ¬ 
ing the coupling of the system to the LTRl in the interaction 
picture reads. 




P ^ 

r s 

I- 


(A14) 


Following the same procedure summarized above, we obtain (m,n - 1,2) is the corresponding lower operator, rid is the 
Eq. Q, where &ki - \ak){ai\ (k,l = 1,2) or amn = \o:m){an\ average phonon numbers of LTRl. The spontaneous decay 
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rates are 

2/7 Vphngklpgmnp 

Jklmn - -■ 

TTC 

with Vphn the phonon volume. The indicator p is defined 
as p - (/ = 1,2). gki,mnp is the coupling constant. 

Jkimnikl = inn) are the decay rates from the upper levels |ai) 
and \a 2 ) to the lower levels |ai) and 1 ^ 2 ), respectively. The 
cross-couplings jkimn (kl + mn) would lead to interferences. 


As the involved levels (coupled to LTRl) here are more than 
those to HTR or LTR2, there are additional cross-couplings 
that do not exist in the case of HTR or LTR2. We obtain two 
types of cross-couplings in the case of LTRl. One is the same 
with that of HTR or LTR2, including 711,21 and 712 , 22 - The 
other type includes 711 , 12 , 711 , 22 , 712 , 21 , 721 , 22 - The first type 
describes transitions to the same lower state, i.e., |ai) for 711,21 
and \a 2 ) for 712 , 22 - But the second type is different. 
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